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Abstract 



iLj ■ We present a detailed analysis of the Connecting Nearest Neigh- 

' . | bors (CNN) model by Vazquez. We show that the degree distribution 

follows a power law, but the scaling exponent can vary with the pa- 
rameter setting. Moreover, the correspondence of the growing version 
of the Connecting Nearest Neighbors (GCNN) model to the particu- 
lar random walk model (PRW model) and recursive search model (RS 
model) is established. 



1 Introduction 

It is well known, that the structure of growing networks is influenced by the 
dynamical processes of their creation such as node addition and edge linking 
[H [21 HI [5] . To show this explicitly, several models of dynamically evolving 
networks have been proposed O HI [6]. 

Many real networks (such as the Internet, for example) expand in time. 
The number of nodes added to the system far exceeds the number of vanish- 
ing nodes. Such networks are well described by the growing models [3]. The 
dynamics of the network growth is captured by the integro - differential dy- 
namical equations. Some growing network models are analytically solvable, 
but often we have to rely on numerical simulations or approximate solutions 
[H HI [5] . Numerical simulations of the growing network usually include 
the following processes: 

• at each time step a new node is added to the network, while bringing 
m new links 

• nodes are labeled according to the time s in which they joined the 
network 

• each link attaches itself to the older nodes with certain probability, for 
example preferentially. 

Barabasi and Albert [3j showed, both numerically and analytically, that 
the preferential node attachment leads to scale free network structure. Pref- 
erential node attachment means that the linking probability of the edge from 
the newly added node to an existing node is proportional to the receiving 
node degree. Ignoring edge direction, the node degree simply counts the 
number of edges incident with that node in the network. The notion of 
"scale free-ness" signifies that there is no dominant scale defining the net- 
work, except of its size. Scale free networks have complex structure with 
degree fluctuations of all sizes, as manifested in the power law degree distri- 
bution 

P(k)ock^, (1) 

where P(k) denotes the probability of choosingj a node with the degree k 
and 7 > is a scaling exponent depending on the network dynamics. It has 
been shown analytically, that for the Barabasi - Albert model ^ba = 3.0 



1 w.r.t. uniform probability over nodes in the network 



[3]. If the node linking differs from the 'pure' preferential one, 7 differs 
from ^iba [DHL or even the network is no longer scale free. For example, if 
the linking probability is random, degree distribution looses its power law 
character [U [5] . 

There are several real examples of scale free networks created by the self 
organizing processes. Some of them, such as the Internet at the autonomous 
system level, social network or language network [2 [10], [TT] have even a 
scale free hierarchical structure. Hierarchy in growing networks has been 
studied by Ravasz and Barabasi p2] and also by us [2J. It has been shown, 
both numerically and analytically, that the hierarchical node organization 
is reflected in the power law distribution of average clustering coefficients 
C(k) of nodes with the degree k 

C{k) oc k- 5 . (2) 

Here 5 > is a scaling exponent. Scale free hierarchical networks have power 
law degree and clustering coefficient distributions (pQ) and ([2]), respectively. 

Pure preferential linking leads to scale free, but not to hierarchical net- 
work structure. It has been shown recently |12} 12]. that both eqs. ([!]) and 
(|2|) hold, if the process of network growth involves some local dynamics - the 
first of the new m links brought by a new node s to the system is attached 
randomly with linking probability 7 (where t ~ N, the number of nodes), 
and the other m — 1 links attach to the neighbors of the first choice node. 
The general scale free hierarchical model (SFH model) has been studied nu- 
merically [21 [12]. For its simplified version with m = 2 (SSFH model) we 
succeeded in finding scaling exponents 7 and 5 analytically [2]. 

A different framework for studying complex networks has been presented 
by Vazquez in [1] - the network structure is discovered through random walk 
by several surfers. Such a random walk model can be transformed to our 
SSFH model [2]. Another model by Vazquez, namely the connecting near- 
est neighbors (CNN) model (section [2]) , has been solved analytically using 
the idea of potential edges and potential degree of the node pp. Differen- 
tial equations, describing the model dynamics, were treated using several 
approximations . 

In this paper we present a more accurate analysis of the the CNN model. 
We show, how the asymptotic degree distribution changes when the line of 
transcritical bifurcation in the parameter space is crossed. We relate our 
more detailed analytical solution to the original solution by Vazquez pQ. We 
also show, that the CNN model corresponds (for certain parameter values) 
to the other known models described in pQ. 



The paper has the following organization: In the second chapter we 
present the CNN model. Chapter 3 is devoted to the qualitative analysis 
of the nonlinear system of differential equations describing the dynamics of 
the CNN network. In the fourth chapter the correspondence of the growing 
version of the CNN model (GCNN model [I]) to the particullar random walk 
model (PRW model [I]) and the recursive search model (RS model [I]) is 
discussed. Section 5 concludes the paper by summarizing the main findings 
of this study. 

2 The CNN model 

The CNN model has been inspired by social networks [lj. It is assumed 
that in such networks two sites with a common neighbor are connected 
with greater probability then two randomly chosen sites. This reflects the 
fact, that in the social network it is more probable that two people (nodes) 
know themselves (are connected) if they have a common friend (common 
neighbor). An analytical understanding of the CNN model has been achieved 
using the notions of potential edges and potential degree. 

The dynamics of the CNN network is defined by the transition rates 
of link states from node s (each possible link from s to other nodes in the 
network can be in three possible states - disconnected (d) , potential edge (p) , 
or an edge (e)). For each node s the transition rates u x ^ y {s) x,y G {d,p,e} 
are defined per link. A potential edge is defined as follows: two nodes are 
connected by a potential edge, if they are not connected by an edge and 
they have at least one common neighbor. For each node s in the growing 
CNN network a degree and a potential degree can therefore be defined. 

Consider a network with N nodes. Denote the (real) degree and potential 
degree of a node s by k(s,N) and k*(s,N), respectively. The dynamics of 
the network is described by [Tj, : 

dk(s,N) 



dN 

dk*{s,N) 



u d ^ e k(s,N) + v p ^ e k*(s,N) - (v e -+ d + v e ^ p )k{s,N) (3) 
Vd-+pK s i N ) +v e ^pk{s,N) - {v p ^ d + v p ^e)k* (s, N) (4) 



dN 

k(s,N) = N-k(s,N)-k*(s,N), (5) 

where fd^e, ^->e depend, as will be shown later, only on the network size 
and Ud^p(s), Vp^d{s) are given as 

Vd->p(s) = U d -> e k(s), Vp^d{s) = Ve^dk(s). 



The last two equations hold since if, for example, a node i' is connected to 
another node i by an edge, potential edges from i' are created to all of the 
neighbors of the node i. 

To simplify eqs. dSHS]) several assumptions reflecting to properties of real 
social networks were introduced [T]: 

1. All processes leading to edge deletion are neglected: f e _> p = zV-><2 = 0. 

2. The transition from a potential edge to an edge has a higher probability 
of occurrence then the transition from being disconnected to an edge 



Mi Mo_ 

N ' Vd ^ e ~ N 2 ' 



where \i\ > and /xo > are constants (for details see [T]). 
3. Terms of order 1/N 2 are omitted. 

Under these conditions the system (J3H5J) turns into 



with solution (&o and /cq are positive constants) [1] 



fc(a,J\0 = *o^-J (8) 

fc*( S ,JV) = fc*(-J , (9) 

where 



o-tO-W^)- (10) 

The degree distribution of the system is therefore power law ([TJ with the 
exponent 7 = 1 + 4 [3 [6] . 



3 Analysis of the CNN model 

In what follows we suggest a more accurate analysis of the original system ([3l- 
[5]) describing the dynamics of the CNN network. From the three simplifying 
assumptions described above and employed in [lj we use only the first two. 
We thus start with the system of differential equations: 



?!%P- = ^(N-k( S ,N)-k*(s,N)) + ^k*(s,N) (11) 

dk * ( d S N N) = ^k( S ,N)(N-k(s,N)-k*(s,N))-^k*(s,N).(12) 

Unlike in [T], our aim here is a qualitative analysis of this nonlinear system. 
We would like to describe asymptotic behavior of solutions of f)l lti!2|) in the 
space of parameters (jlq, fi\. The goal of the following series of substitutions 
is to transform the system (jlltil2p to a more convenient form, which can be 
approximated by an autonomous nonlinear system. 

In the continuum approach [6] N is regarded a continuous variable 
and s is a node label. First set N = e T and substitute 



x(s,t) = k(s,T) + k*(s,T) (13) 

y(s,T) = k(s,r) + l, (14) 

to the system (fTTHT2|) . We obtain 

x = /i (l - xeT T )y (15) 

y = /i (l -xeT T ) + m(x-y + 1), (16) 

where x(s, r) and y(s, r) are abreviated as x and y, respectively and we 
write x' and y' for X q 1 1 t and y g' T , respectively. 
Using another substitution 



u = l - xe - T (17) 

v = ye-\ (18) 

together with abbreviations u(s,t) = u and v(s,t) = v, we again get a new 
system of differential equations: 



u = 1 — u — fiQUv (19) 

v' = [j,i — \x\u — (jl\v — v + iJ>oue~ T + Hie~ T , (20) 

with obvious interpretation of v! and if. It is natural to assume that the 
function u(s,t) (fTTl) is boundeco Then for r — > oo, the functions Ho ue ~ T \ 
Hie~ T tend to zero. We approximate the system (fT9H20|) by an autonomous 
one 

v! = 1 — u — fiouv (21) 

v' = Hi — fiiu — Hi v ~ v - (22) 

The correctness of the approximation (|21H22p can also be shown if we intro- 
duce a new function z = e~ T and rewrite the system (|19H20p into a three 
dimensional one 

u = 1 — u — hquv (23) 

v' = Hi ~ f J 'i u ~ A*i u — v + IMsUZ + Hi z - (24) 

z 1 = -z. (25) 

Now, we can see from (|23H25p that all of the important dynamics occurs in 
the (u, v) plane. 

The system (|2TH22|) possesses two fixed points 

A = [1,0] (26) 

and 

B = [ /il + 1 t2tl -//1 ~ 1 1 (27) 

I MoA*i ' /"o/"i + Mo J 

Having identified the fixed points, the next step is to characterize their 
stability type. To that end we linearize (|21H22|) . The linearization matrix 
M of {SUisH reads 

M= (-l-^v -hou \ (2g) 

\ -Hi -Hi ~ 1 J 



2 As the network size N = e T grows, the degrees (both actual and potential) of any 
node s < N cannot grow faster than c • iV for some positive constant c. 



and at the fixed point A ([26]) we get 

M A =( ~ l '"' , I. (29) 




Analogously at the fixed point B (127p we get 

w+i mi (30) 

-/il -fil-1 J 

The types of fixed points A, i? depend on eigenvalues of Ma, Mb- General 
form of eigenvalues of the matrix M (|28p is 

where Tr(M) and Det(M) denote trace and determinant of M, respectively. 
Clearly 



^^ = -l-f, Det(M A ) = l + , 1 -^ 1 



and 



Tr(M B ) fipfn 

r =-l-yUi — -, Det(M B ) = -1 - fil + fiQfil- 

I fJLi + 1 

Both traces Tr(M^) and Tr(Ms) are negative for each /j,q > 0, /xi > 
(which, of course, is fulfilled for the CNN model), and the determinants have 
opposite signs, Det{MA) = —Det(Mg). It therefore follows from (|3ip that 
the sign of the expression 1 + /xi — ^o^i is decisive for the stability type 
of A and B. In the case that 1 + m — /^oW > 0, the eigenvalues of M^ 
are both negative and A is an asymptotically stable fixed point called sink. 
The eigenvalues of Mb have opposite signs and B is an unstable fixed point 
(saddle). In the complementary case, 1 + fi\ — ^o/^i < 0, A and B change 
their stability types, now A is a saddle and B is a stable fixed point (sink). 
The critical case occurs when 1 + /ii — /io^i = 0. In that case the fixed 
points A, B coincide. 

This type of stability change is called a transcritical bifurcation (e.g. 
[8, 9j). In the CNN model the bifurcation occurs when the parameters 
(/xc^i) cross the line 

Mo = 1 + — (32) 

Mi 



A saddle, B sink 



A sink, B saddle 



Figure 1: Bifurcation line obtained from (|32D in the parameter space of 
the CNN model delimiting regions of different stability types of fixed points 
of (jl!H12p . Dotted horizontal line indicates an asymptotic behavior when 

\x\ — > oo. 



in the parameter space (see Fig. [[]). 

Let us now consider both cases in greater detail. In the first case, 1 + 
fj-i — fJ-ofJ-i > 0, A is a stable sink and the eigenvalues of Ma are 



Ai, 2 = -l - y ± yjw + y- ( 33 ) 

The corresponding eigenvectors of Ma for eigenvalues Ai, A2 read 



Wl 



''2 



-Mo 



/'I 

2 


+ V M0M1 + — 
-Mo 


A'l 

2 


- V M0M1 + — 



(34) 
(35) 



The general solution of the linearization which approximates solutions of the 
system ()21H22p and therefore also solutions of ()19H20p near the stable fixed 
point A is 



u 

v 



' I + ciUie AlT + c 2 v 2 e X2T ■ 







(36) 



where c\, c 2 are constants. Returning to the original variables, neglecting 
e A 2 T t erms anc j U si n g initial condition k(s,s) = 1 we get the final approxi- 
mation 



k(s,N) 
k*(s,N) 
where the exponent j3 is given by 



'-% 



N\P 



1 



"$->)&+*• 



(37) 
(38) 



-1 + 1 1 + 



4^o 
Mi 



(39) 



Equation (|37p allows us to calculate the average degree as 



k(N) 



1 N 1 



s=l 



N 



s<r 



i 



(40) 



Approximating the sum in (|40p by definite integral, the asymptotic (as N 
oo) average degree reads 



1 + 13 
1-/3' 



(41) 



There is another way of obtaining the asymptotic average degree k. From 

(EHD, 



k (s,s) = 2— -1 = WH . 

P V Mi 



(42) 



This result has an interesting interpretation. If a node comes and links itself 
to a node s' by one edge, potential edges are immediately created to all of 
the neighbors of s'. Therefore k*(s,s) can be interpreted as an asymptotic 
average degree. The difference between (|4"Tj) and (|4"2"|) is due to the neglected 
terms in (f37l) . but both calculations lead to the finite average degree k. 



10 



In the other case, when 1+/j,i — /UoMl < O7 fixed point B is asymptotically 
stable and the eigenvalues of Mb are 



l + /ii // /ii //1 + Mi . M0M1 \ 2 , . 1 

(43) 
Clearly A 4 < A 3 < 0. 

The general solution approximating solutions of the system (|19H20p near the 
stable fixed point B is therefore 

f \ / Ml + l \ 



■y / I Wf 



^-1 + C3V 3 e MT + c 4 v 4 e MT , (44) 

MOM1+M0 / 

where C3, C4 are constants and v% = (vn,Vi2) , i G {3,4} are eigenvectors of 
M B - Let us denote a = MoW ~^ 1 " 1 . 

Backward substitutions iead to the solution in the original variables k, k*: 



k( s , N) = fj,N - c 3 v 32 N 1+X3 - c 4 v 42 N 1+X4 - 1 (45) 

k*(s,N) = —N-c 3 (v 31 -v 3 2)N 1+X:i -c 4 (v 41 -v 42 )N 1+Xi + l(46) 
Mi 

where coefficients c 3 c 4 are s -dependent. 

Since A4 < A3 < 0, the largest exponent of N is 1 and we can approximate 
the solution by 



k(N) = nN (47) 

k*(N) = -^-N (48) 

Mi 

For N — > 00, 



lim H *' N \=^ (49) 

and /Lti is thus an asymptotic ratio of the node's degree and its potential 
degree. 

It is possible to find a neater approximation. As A4 < —1, neglecting 
only terms with A^ 1+A4 in (|45p we obtain 
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k(s,N) = fiN - c 3 v 32 N 1+X:i - 1 (50) 

k*(s,N) = ^N-c 3 (v 31 -v 32 )N 1+x *+l (51) 

Ml 

which together with the initial condition k(s, s) = 1 leads to the final result 



fc( S ,iV) = ^^l-^^^Ar + 2^ 



1 (52) 



fc*(*,JV) = ^('l+(^) A3 )iV-2^) 1+A3 - C3U 3iiV 1 + A 3 + 1( 53 ) 



where 

U3i = 1 + — • (54) 

Ml 

Again, from (|52p we can derive the average degree: 

1 N 1 N { /M\ X A //V\ 1+A » 

^.-gK^J.-g^!^-) )n + 2(-) -1.(55) 

Approximation of the sum in (|55|) by definite integral and taking the limit 
N — > oo leads to 

k = co. (56) 

As before, the same asymptotic result can be obtained from 

k* (s, s) = —2s - 1 - c 3 v 31 s 1+Xi (57) 

Mi 

and 

Urn k*(s,s) = k* = oo. (58) 

s— >oo 

The asymptotic average degree tends to infinity when the fixed point B is 
stable. 

To conclude our analysis we can state that the CNN model ([3][5]) exhibits 
two different types of behavior, deliminated by the bifurcation line (|32p in 
the parameter space (//oiMi)- 

12 



For fiQ < 1 + — the degree k(s, N) with increasing N grows sublinearly 
with the exponent /3 given by (j39|) , where < /3 < 1 . Due to approximations 
(|4ip . (|42p the asymptotic average degree is in this case constant. The degree 
distribution is a power law with exponent 7 = 1 + 4 [5j [6]. 

For /Uq > 1 + 7p, the growth of the degree k(s,N) is linear (exponent 



13 = 1) with increasing iV. Approximations (|56p . (|58p imply that now the 
asymptotic average degree tends to infinity. The degree distribution is a 
power law with the exponent 7 = 2. 

Table 1: 7 exponent of the degree distribution in a parameter space region 
where the fixed point A is a sink and B is a saddle. /j,q = 1, 74 is calculated 
with a help of (jlOl [39]) and 7„ is a numerical value. 



j»l 7t 7n 

0.25 3.56 3.49 

0.5 3.0 3.03 

1.0 2.62 2.66 



We supported our analytical results by the numerical simulations. We 
simulated the model for parameters /io, fi± in both regions of the parameter 
space (Fig. [[]), using the equations (fTTl [T2"P . Our network has 100000 nodes. 
In the first region, where the Vazquez solution holds, we simulated the model 
for no = 1 and for the three different values of [X\. The numerical results 
are in a good accordance with the theoreticaly predicted values (Tabfl]). 

Above the bifurcation line the numerical simulations does not converge 
very well for the parameter values no, [i\ close to the bifurcation line. We 
belive, this is partly due to the fact, that the second terms of the equations 
(|45p and (|46|) influence the convergence. The results for several parameter 
values are in the (Tabl2]). (Fig. [2]) shows the convergence of the 7 exponent 
for /xq = 2, 3 and growing [i\. 

4 Correspondence of the CNN model to other 
models 

Having obtained analytical solutions for the CNN model, we would like 
to establish whether some known models are, at least for some parameter 
settings, identical to the CNN model. We have tested a growing variant of 
the CNN model (the GCNN model) and the recursive search (RS) model, 
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Table 2: 7 exponent of the degree distribution in a parameter space region 
where the fixed point A is a saddle and B is a sink, jt is a teoretical value 
of the scaling exponent, and 7 n is a numerical value. 



fJ-o j»l It In 

2 6 2.0 2.21 

2 8 2.0 2.08 

3 9 2.0 2.19 
3 12 2.0 2.06 



which is derived from the PRW model. All of them are described in pQ. 

4.1 The GCNN model 

In pQ Vazquez proposed a growing variant of the CNN model, the GCNN 
model. Fix some u G (0, 1). The GCNN network evolves by following 
iteratively performed rules: 

1. With probability 1 — u a new vertex is introduced and a new edge 
from the new vertex to a randomly selected network node is created. 

2. With the probability u one potential edge (selected at random) is 
converted to an edge. 

The author expects that in the regime of large N, the CNN and GCNN 
models will behave similarly. It has been stated that the evolution rules are 
consistent with setting /iq = 1 and fii = j^ in ((6]-[7|): 



«K.,.v) _ 1 / u s 1 



dN N \l-uj N 

3 k-(s,N) H .,N) I u ^ 



dN N \l-uj N 

How to interpret this? The number of nodes N is now proportional to the 
time of the network development. If the time is measured by coming nodes, 
it is necessary to rescale k* by the factor j^ which is, in fact, the ratio of 
transformed potential edges to one node. For our next analysis it is quite 
correct to deal with the original set of equations ([6HZ1 l59"H60|) , for reasons 
which are to be made clear later. 
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3.2 



^-2.6 




Figure 2: Convergence of the 7 exponent to the predicted value 7 = 2.0 
in the region above the line of the transcritical bifurcation. /j,q = 3 for the 
upper curve (x) and /Uo = 2 for the lower curve (+). All simulations were 
done for 100000 nodes. 



Having the parameter settings /j,q = 1 and /ii = j^— and taking into 
account that 



7 = 1 + 



ir 



(61) 



the scaling exponents (3 and 7 of the GCNN model are 



(3 



a 



2(1 -u) 



I + 1/I + 4- 



1-u 



(62) 



and 



2(1 -it)/ / 1-u 
7 («) = 1 + -^ '- [ -1 + a/1 + 4 



11 



it 



(63) 
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The limiting cases, analyzed in pQ are 

7(0) = 00 (64) 



and 



7 (1) = 2. (65) 



Viewing these results in the light of our solution, it is clear that the condition 
fj-o = 1 ensures, we are in the region of the parameter space, where the fixed 
point A is an asymptotically stable sink (Fig. [1]). Exponents /3 and 7 are 
therefore given by (|62p and (|63p . This is also why it is correct to use the 
original equation set ([H \59\i , which gives a good result in this part of the 
parameter space. From ([63]) we have that 7(0) = 00. 

There is an alternative argument: If u = 0, only the node addition rule 
works, each new node creating a new edge to a randomly selected node 
in the network. Such network, growing by the random node addition, has 
been analyzed in [5] E] • In such networks the degree distribution decreases 
exponentially, which corresponds to 7(0) = 00. 

To analyze the second limiting case, we rewrite (|32p . 

Mi = ~^—r- (66) 

Mo - 1 

Comparing (|66p with [i\ = t^— , it is clear, that for — = and u = 1 both 
expressions are identical. That means, that for the limiting case u = 1, 
we are assymptotcally on the bifurcation line (|32j) in the parameter space, 
where the solution changes to that described by the equations (|47p . Here 
the scaling exponents are j3 = 1 and 7 = 2. 

Different scaling exponents in the two cases 7(0) and 7(1) (|63[) are there- 
fore a natural consequence of the fact, that the solution lies in different 
regions of the parameter space divided by the line of the transcritical bifur- 
cation (l32ll. 



4.2 The GCNN model: u=\ case 

We get another interesting result when (3 = | ([62]) . Using ([61]) leads to 
7 = 3. This value of the 7 exponent is typical for several well known complex 
network models, for example the Barabasi Albert model (BA model) [3] or 
the scale free hierarchical model (SFH model) [2]. In [2] we established a 
simple process of network growth that leads to scale free hierarchical network 
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structure in which both degree and clustering coefficient distributions are 
power law (HJ [2]) . 
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Figure 3: Log-log plot of the degree and potential degree developement with 
the size of the GCNN network for two different nodes (s = 5000,10000). 
k(s,N) (k*(s,N)) is a degree (potential degree) of the node s and network 
size N. The line has a slope \ (|7Tj) . 



In the GCNN model either new vertex with new edge is added to the 
system with probability 1 — u, or, with probability u, some potential edge 
is chosen and turned into an edge connecting two nodes with a common 
third node to which the two nodes are connected. A triangle is thus created. 
Node addition causes creation of several new potential edges in the system. 
Let us assume that 



dk*(s,N) dk{s,N) 



(67) 



dN ON 

holds for the GCNN model for some value of u. Using this assumption, we 
get from ([59]) 
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«•(..*)_!/ - U„ n (68) 



<9iV Af \1-u/jV 

with the solution 

k*(s,N)=[—) c-- -. (69) 

V s ) u 

where the exponent f3 (u) = j^— and c is a constant. Comparing this value 
of /3 with ([62]) we get that u = |. That means, that /3 ( |) = | is a result 
of the fact, that, as the network grows, both degree and potential degree of 
nodes change the same way. 

Summing ([59]) and ([60]) and using the assumption ([6?]) we get 

2 ™4< i+ ^> po) 

with solution (under initial condition k(s,s) = m) (Figj3]) 



/c(s,iV) = (m+l)(— J -1 (71) 

i 



fc*(s,JV) = (m + l) f — J -2. 



To conclude this analysis, we can state that if the special condition 
— 9 ^y ' = — ^ ' holds in the GCNN model, then u = | and the scalings 
exponents f3 = g and 7 = 3.0. It would be interesting to investigate, whether 
the GCNN model can be transformed to some other known models having 
degree distribution governed by the exponent 7 = 3.0. We belive, that the 
scale free hierarchical (SFH) model [2] and the GCNN model are similar, 
but we leave this to future work. 

4.3 Correspondence of the GCNN model to the RS model 

Recursive search model (RS model) has been investigated by Vazquez in 
PQ. It is a variant of his particular random walk model (PRW). General 
random walk model simulates how the structure of a network is discovered 
by surfers pQ. The model has been inspired by the manner people obtain 
information by surfing on the world wide web. Here the web pages are nodes, 
and hyperlinks are edges. 
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How does the surfer move? There are two possibilities. Either to jump 
on a page selected randomly, and then follow a hyperlink, or to jump on 
another randomly selected page. PRW model is a simpler version of this 
random walk model. The dynamics of the PRW model consists of adding 
and walking. The process starts with a single node. Adding means, that 
the new vertex is created and connected to the one of the existing vertices 
by an edge. Walking is defined as follows: if an edge is linked to a vertex 
in the network, the edge is also created to one of its nearest neighbors with 
probability q e . When no edge is created, adding rule is applied. 

In the RS model the walking rule is slightly changed [I]. Here walking 
means that if an edge is created to a vertex in the network, then with 
probability q e an edge is also created to all of its nearest neighbors. When 
no edge is created the adding rule is applied with the probability 1 — q e . 

Vazquez pQ solved the RS model analytically for q e = and q e = 1. 
In the first case (q e = 0) only adding is present. No link is given to the 
neighbors of the chosen vertex. We believe that this is the same situation as 
in the GCNN model with u = 0. Here adding creates new potential edges, 
but none of them is turned to an edge. It has been shown, that the degree 
distribution in the recursive search model decreases exponentially, which 
corresponds to the 7(0) = 00 (j63j) case in the GCNN model. And, again, 
we can say, that the network dynamics is the same as that of the network 
growing by the random node attachment [U [3] . 

In the second limiting case q e is close to 1 and adding in the RS model 
is almost not present. But if a node is eventually added, link is created to a 
randomly chosen node in the network and "potential edges" among the new 
node and all of the neighbors of chosen node are created and immediately 
converted to edges. This corresponds to the GCNN model with u = 1. Each 
time step a potential edge is transformed to an edge, with almost no adding, 
all potential edges are finally transformed into edges. Degree distribution 
for both models is a power law with scaling exponent 7 = 2. 

5 Conclusion 

We present a detailed analysis of the CNN network model [lj. The parame- 
ter space is divided into two regions in which the CNN model has different 
behaviour (FigQ]). The degree distribution in both regions is a power law, 
but the scaling exponent 7 is different in different parts of the parameter 
space. We supported our analytical studies by the numerical simulations. 
While the results belove the bifurcatin line (Fig(T]) closely follow the pre- 
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dieted values (Table. 1), above the line the convergence to the predicted 
value 7 = 2.0 is worse in the proximity of the bifutcation line (Table. 2), 
(Figj2]). We have nevertheless shown, that the 7 exponent converges to the 
theoretically predicted value. 

Our theoretical analysis enabled us to understand the dynamics of the 
growing variant of the CNN model, the GCNN model, in two limiting cases 
of u = and u = 1. The two cases belong to the different regions of the 
CNN parameter space. 

We analysed the GCNN model with u = 3, in which case the scaling 
exponent 7 of the degree distribution is 7 = 3.0. We have found, that this 
results from the fact, that both the degree and the potential degree of nodes 
in the GCCN network changes similarly. We speculate that the same holds 
for the SFH model and that the two models are related when u = 3. 

Moreover, we have shown that the RS model and GCNN model in the 
limiting cases u = 0, 1 and q e = 0, 1 have similar dynamics. 
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